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We present measurements of direct photon pair production cross sections using 8.5 fb _1 of data 
collected with the DO detector at the Fermilab Tevatron pp collider. The results are presented 
as differential distributions of the photon pair invariant mass dcr/dM 77 , pair transverse momen- 
tum dcr/dpj? , azimuthal angle between the photons dcr/d A</> 77 , and polar scattering angle in the 
Collins-Soper frame da/d\ cos#*|. Measurements are performed for isolated photons with trans- 
verse momenta > 18 (17) GeV for the leading (next-to-leading) photon in pr, pseudorapidities 
|?7 7 | < 0.9, and a separation in n — (j> space A7Z~ n > 0.4. We present comparisons with the predic- 
tions from Monte Carlo event generators diphox and RESBOS implementing QCD calculations at 
next-to-leading order, 27NNLO at next-to-next-to-leading order, and SHERPA using matrix elements 
with higher-order real emissions matched to parton shower. 
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Precise knowledge of the direct diphoton (DDP) pro- 
duction differential cross section is a cornerstone of the 
search for the standard model (SM) Higgs boson by ex- 
periments at the Large Hadron Collider [l], 0] and the 
Tevatron 0, Q|. The term "direct" means that these 
photons do not result from meson, for example, 7r ,7y,u;, 
or Kg decays. DDP production is also a significant 
background in searches for Kaluza-Klein [|| or Randall- 
Sundrum [|| gravitons decaying into two photons, as well 
as other new phenomena processes, such as decays of 
heavy resonances Q or cascade decays of supersymmet- 
ric particles Q. For these searches, DDP production is 
an irreducible background, and it is crucial to have a de- 
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tailed understanding of the distributions of key kinematic 
variables [§]. 

In addition to investigating physics beyond the SM, 
DDP production processes are important for studying 
quantum chromodynamics (QCD) and measuring parton 
distribution functions (PDFs). DDP production cross 



sections have been examined at fixed-target jlOt lllj and 
collider experiments (l2l - fl7j . DDP events at the Tevatron 
pp collider are produced predominantly through quark- 
antiquark annihilation qq — > 77 and gluon-gluon fusion 
(.9.9 "~ * 77) v i a a quark-loop diagram. The matrix element 
(ME) for the latter process is suppressed by o? s relative 
to qq annihilation, but its total production rate at low 
77 invariant mass (Af 77 ) and intermediate 77 transverse 
momentum (jhf 1 ) is quite significant due to the relatively 
large values of the gluon PDFs in that kinematic region. 
By the same argument, gluon-gluon fusion becomes even 
more important at the LHC [l8j. DDP events may also 
originate from processes such as qg — > 97, 99 — > .97, and 
.9.9 — ^ 99S where a photon with large transverse momen- 
tum is radiated from the final state parton. These pro- 
cesses, being nearly collinear, require the introduction of 
a fragmentation function in perturbativc QCD (pQCD) 
calculations 18 1. Photon isolation requirements reduce 
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the contribution of such fragmentation events. How- 
ever, their contribution may be still quite large at low 77 
azimuthal angle difference (A$> 77 ) and for intermediate 



77 



19( , which is the DDP transverse momentum. 
In this Letter, we present measurements of differential 
cross sections of DDP production using the dataset col- 
lected at the Fcrmilab Tcvatron DO experiment between 
June 2006 and September 2011. The dataset corresponds 
to an integrated luminosity of 8.5±0.5 fb -1 [35j . 

Measurements are performed as functions of M 77 , 



A0 77 , and |cos0*|, the absolute value of the co- 



77 
Pt . 

sine of the polar scattering angle of the diphoton sys- 
tem in the Collins-Sopcr frame [23]. Here we approxi- 
mate I cos 0*| by I tanh[(?7i — T72 ) / 2] | , where 771,2 are the 
pseudo-rapidities 2l| of the leading and next-to-leading 
photons ranked by px- These four variables emphasize 
different phenomena in the diphoton production mecha- 
nism. M™ usually serves as a probe for new phenomena 
searches [lj, [4, l5j-|7| and PDFs. The p T 7 and A</> 77 shapes 
are mostly sensitive to the initial state gluon radiation 
and fragmentation effects. The | cos 8*\ angle is sensitive 
to PDFs and spin correlations in the final state. In con- 
trast with the previous DO measurement [l4| , in this anal- 
ysis we do not impose explicit minimum requirements on 
M 77 or A</> 77 , nor do we require that M 77 > p^ 7 , mak- 
ing the measurements more universal. By separating the 
data into two subsets, with A</> 77 > n/2 and A<^> 77 < tt/2, 
we isolate regions with smaller and larger expected rela- 
tive contributions from the fragmentation processes. 

We compare our results with the theoretical predic- 
tions generated using the DIPHOX [lU, RESBOS (9I I22I [23| . 
27NNLO 24 1 and SHERPA [25| event generators. The gen- 
eral multipurpose generator approach is to employ inter- 
leaved QCD and quantum electrodynamics (QED) par- 
ton shower (PS) to describe initial and final state radi- 
ation. The SHERPA Monte Carlo (MC) event generator 
improves this technique by including higher-order real- 
emission matrix elements pGj . Matching between par- 
tons coming from real emissions in the ME and jets from 
PS is done at some (hardness) scale Qc U t defined follow- 
ing the prescriptions given in Ref. [26[. We use events 
generated with all MEs with two photons and up to two 
hard partons. However, the ME for gluon-gluon scat- 
tering gg —> 77 in SHERPA does not have real parton 
emissions. As is shown in Ref. 26j, SHERPA provides 
a good description of the fragmentation function mea- 
sured at LEP at high fraction of the jet energy carried 
by the photon, corresponding to tight photon isolation 
cuts. The loop corrections matching the higher order 
MEs are missing in SHERPA, which can make predictions 
significantly scale-dependent and may lead to underesti- 
mation of 77 rates. In the SHERPA version used in this 
paper p5| , the inherent next-to- leading-logarithmic effect 
of correlated emissions is invoked in parton-shower sim- 
ulations by appropriately choosing a scale factor for the 
argument of the running strong coupling constant [27i — 



l29(. The DIPHOX and RESBOS packages provide predic- 
tions at next-to-leading order (NLO) in pQCD, with the 
99 ~~ ^ 77 process considered only at the leading order 
approximation in diphox. Also, in diphox, explicit sin- 
gle and double parton-to-photon fragmentation processes 
are included at NLO accuracy, while in RESBOS, rates of 
fragmentation processes are approximated by a function. 
Only in RESBOS are the effects of soft and collinear initial 
state gluon emissions resummed to all orders (23|. The 
rcsummation should be important for a correct descrip- 
tion of the p T 7 distribution close to zero and the A0 77 
distribution close to tt. The 27NNLO generator, which 
appeared recently, exploits the pip subtraction formalism 
[30l | that handles the unphysical infra-red divergences up 
to next-to-next-to- leading order (NNLO). It takes into 
account most diagrams (qq and qg scatterings) at 0{a 2 s ) 
accuracy; however, in the current calculations, there is 
no higher-order correction to the gg — > 77 box diagram 
and no soft gluon rcsummation is applied. Additionally, 
it does not take into account the fragmentation contri- 
butions. 

The DO detector, where the DDP measurements are 
performed, is a general purpose detector described in de- 



tail elsewhere 31-331. The sub-detectors used in this 



analysis to trigger events and reconstruct photons are the 
calorimeter, the central tracking system, and the central 
preshower. The muon detection system is used to com- 
pare data and MC simulation sets of Z — > + 7 
events to obtain data-to-MC scale factors for reconstruc- 
tion efficiency. The central tracking system, used to re- 
construct tracks of charged particles, consists of a silicon 
micro-strip detector (SMT) and a central fiber track de- 
tector (CFT), both embedded in a 2 T solcnoidal mag- 
netic field. The solenoid is surrounded by the central 
preshower (CPS) detector located immediately before the 
inner layer of the electromagnetic calorimeter. The CPS 
consists of approximately one radiation length of lead ab- 
sorber surrounded by three layers of scintillating strips. 
The calorimeter is composed of three sections: a central 
section covering the range of pseudo-rapidities |f?dct| < 
1.1 [2l| and two end calorimeters (EC) with coverage ex- 
tending to |r?det| ~ 4.2, with all three housed in separate 
cryostats. The electromagnetic (EM) calorimeter is com- 
posed of four layers of A?7dct x A0dct = 0.1 x 0.1 cells, 
with the exception of layer three (EM3) with 0.05 x 0.05 
granularity. The calorimeter resolution for measurements 
of the electron/photon energy at 50 GeV is about 3.6%. 
The energy response of the calorimeter to photons is 
calibrated using electrons from Z boson decays. Since 
electrons and photons shower differently in matter, addi- 
tional corrections as a function of 77 are derived using a 
detailed GEANT-based [HI simulation of the DO detector 
response. These corrections are the largest, (2.0-2.5)%, 
at low photon energies (« 20 GeV). Events satisfying 
the following trigger requirements are recorded: at least 
two clusters of energy in the EM calorimeter with a loose 
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shower shape requirement and a range of pr thresholds 
between 15 GeV and 25 GeV. Luminosity is measured 
using plastic scintillator arrays placed in front of the EC 
cryostats. 

Events are selected with at least two photon candidates 
with transverse momentum pr > 18 (17) GeV for the 
leading (next-to-leading) candidate and pseudorapidity 
\t)\ < 0.9. We require a slight difference between the pr 
cutoffs for the two photons to avoid a divergent kinematic 
region of the NLO calculations 18| . The trigger is more 
than 90% efficient for these selections. 

At high instantaneous luminosities there is more than 
one pp interaction per beam crossing. The photon pt is 
computed with respect to the reconstructed pp interac- 
tion vertex with the highest number of associated tracks, 
called the event vertex [l4| . The event vertex is required 
to be reconstructed within 60 cm of the center of the 
detector along the beam axis (z), and satisfies this re- 
quirement in 98% of events. We correct for effects of 
selecting an incorrect vertex (in about 35% of events) us- 
ing Z —> e + e~ data events, where we remove tracks cor- 
responding to the electron and positron to model DPP 
production. 

Photon candidates arc formed from calorimeter towers 
m a cone of radius TZ = ^(Arj) 2 + (Ac/)) 2 =0.4 around 



a seed tower [3l|. A stable cone is found iteratively, and 
the final cluster energy is recalculated from the inner core 
within TZ = 0.2. The photon candidates are required to: 
(i) have > 97% of the cluster energy deposited in the 
EM calorimeter layers; (ii) be isolated in the calorime- 
ter according to [£ t ot(0.4)-£EM(0.2)]/£ E M(0.2) < 0.07, 
where E tot (TZ) [E EM (TZ)] is the total [EM only] energy in 
a cone of radius TZ; (iii) have the scalar sum of pt's of 
all tracks originating from the event vertex in an annu- 
lus of 0.05 < TZ < 0.4 around the EM cluster less than 
1.5 GeV; and (iv) have an energy- weighted EM shower 
width consistent with that expected for an electromag- 
netic shower. To suppress electron misidentification as 
photons, the EM clusters are required to have no spatial 
match to a charged particle track or any tracker hit con- 
figuration consistent with an electron. The two photon 
EM clusters are required to be separated by A7\L 77 > 0.4. 

Additional group of variables exploiting the differ- 
ences between the photon-initiated and jet activity in 
the EM calorimeter and the tracker are combined into 
an artificial neural network (NN) to further reject jet 
background [3^|. In these background events, photons 
are mainly produced from decays of energetic ir° and_ 
mesons. The NN is trained on 7 and jet pythia 
MC samples. The generated MC events are processed 
through a GEANT-based simulation of the DO detector. 
Simulated events are overlaid with data events from ran- 
dom pp crossings to properly model the effects of multiple 
pp interactions and detector noise in data. Care is taken 
to ensure that the luminosity distribution in the overlay 
events is similar to the data used in the analysis. MC 



events are then processed through the same reconstruc- 
tion procedure as the data. MC events are reweighted to 
take into account the trigger efficiency in data, and small 
observed differences in instantaneous luminosity and dis- 
tribution of the z coordinate of the event vertex. Pho- 
ton radiation from charged leptons in Z boson decays 
(Z — > l + l~7, £ = e, fx) is used to validate the NN per- 
formance [38[ . The NN describes the data well and gives 
significant extra discrimination against jets. The pho- 
ton candidates in this analysis are chosen such that their 
NN output requirement retains 98% of photons and re- 
jects 40% of jets beyond the rejection provided by the 
selection described above fl4| . 

We estimate contributions from instrumental 7+jet 
and dijet backgrounds and also the contribution from 
Z boson/Drcll-Yan production events Z/7* — > e + e~ 
(ZDY). In the instrumental backgrounds, one or more 
jets are misidentified as photons from jet-forming par- 
tons that hadronize into isolated neutral meson(s) (tt° or 
rf) giving rise to two or more photons in the final state. 
Electrons in the ZDY background can be misidentified 
as photons due to similarities in the shower shape. The 
contribution from the ZDY events is estimated from MC 
simulation with pythia, normalized to the NNLO cross 
section [3^]. On average, 2% of the electrons survive the 
selection criteria above, mainly due to the inefficiency of 
matching a charged track to an electron. In data this in- 
efficiency is higher than in MC and the ZDY contribution 
is corrected for these differences. 

The 7+jet and dijet instrumental backgrounds are es- 
timated by fitting a two-dimensional (2D) distribution of 
the leading and next-to-leading photon NN outputs with 
templates extracted from DDP SHERPA signal and EM- 
jct pythia MC samples. In the latter, constraints are 
placed at the generator level to increase the statistics of 
jet events fluctuating into EM-like objects |2f|. For the 
7+jet template, the photon candidate is taken from ei- 
ther the 77 sample or from the EM-jet sample, while for 
the dijet template, both candidates are taken from the 
EM-jet sample. Table U shows the numbers of events sur- 
viving the selection in data for different A0 77 regions, as 
well as the number of data events from each of the four 
sources as determined by a fit of the signal and back- 
ground templates to data. The typical DDP purity in 
the selected data events is around 60%. 



TABLE I: The numbers of 77 (Y 77 ), T j+j 7 (Y 7J ), jj (Njj), 
and ZDY (JVzdy) events and their total. The quoted uncer- 
tainties are statistical only and for iV 77 , N 1 j, and Njj are 



from 2D fitting. 



Full A0 77 A<£ 77 < 7r/2 A0 77 > 7t/2 
A/" 77 20255+398 1676+109 18572+370 
N yj 2575+516 317+148 2217+459 
Njj 10992+344 854+96 10185+314 
AW 198+14 2.7+1.7 195+13 
Total 34020 2851 31169 
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The estimated numbers of DDP events in each bin are 
corrected for the geometric and kinematic acceptance 
of the photon, as well as for the photon detection effi- 
ciency. Both acceptance and efficiency are calculated us- 
ing SHERPA MC events. The acceptance is calculated for 
the events satisfying at the particle level [Io[ p T > 18 (17) 
GeV for the leading (next-to-leading) photon, |?7 7 | < 0.9, 
and A1Z 17 > 0.4. The photon is also required to be iso- 
lated by p T ° = p)^ (0.4) -j? T < 2.5 GeV, where pJp*(0.4) 
is the scalar sum of the transverse momenta of the sta- 
ble particles within a cone of radius 1Z = 0.4 centered 
on the photon. The acceptance is driven by selection re- 
quirements in ?7dct (applied to avoid edge effects in the 
calorimeter regions used for the measurement) and 
(to avoid periodic calorimeter module boundaries) [3l|, 
photon rapidity rf and pt, and bin-to-bin migration ef- 
fects due to the finite energy and angular resolution of the 
EM calorimeter. Typically, greater than 80% of events 
at the reconstruction level remain in the same bin as at 
the particle level. Choice of an incorrect event vertex 
leads to a systematic uncertainty on the acceptance, typ- 
ically < 3% for A</> 77 > tt/2 and < 6% for A0 77 < tt/2. 
The systematic uncertainty is estimated by using DDP 
events simulated with SHERPA in which the event vertex 
position is randomized according to its distribution in 
z with respect to the true vertex, and by recalculating 
all relevant variables of the diphoton system. Possible 
mo del- dependent effects are corrected by recalculating 
the acceptance according to the difference between the 
photon px spectra in data and SHERPA MC. The accep- 
tance grows from 45% in the low M 77 region to 80% in 
the high mass region. The systematic uncertainty on the 
acceptance varies within (4-21)%. In the regions domi- 
nated by fragmentation photons, such as low A0 77 and 
intermediate p 7 -, 7 , the acceptance is lower than in the re- 
gions dominated by direct production. The EM clusters 
reconstructed in the acceptance region are required to 
pass the photon identification criteria listed above. Small 
differences between the photon identification efficiencies 
in data and MC are corrected by using control samples of 
electrons from Z boson decays and photons from radia- 
tive Z boson decays [3, [HI . The overall diphoton selec- 
tion efficiency is typically about 50% with variations of 
±5%. The relative systematic uncertainty of the dipho- 
ton selection efficiency is about 4%. 

The differential cross sections da/dM 11 , <kr/dp T 1 , 
da I dA0 77 , and da /d\ cos 9* | arc calculated from the 
number of data events after the subtraction of back- 
ground contributions divided by the event selection ef- 
ficiencies, acceptance, integrated luminosity, and the bin 
width. 

The measured differential cross sections for all consid- 
ered kinematic regions are presented in Tables ITll lVl The 
average value of each variable in a bin was estimated us- 
ing SHERPA MC events. The statistical uncertainty <5 s t a t 
is caused by finite MC statistics used for the efficiency 



TABLE II: The measured differential cross sections in bins of 
M 77 and . The columns <5 s tat, S syst , and 5 to t represent the 
statistical, systematic, and total uncertainties, respectively. 



M 77 


<M 77 ) 


do 1 dM-f-f 


5stat 






<5tot 


(GeV) 


(GeV) 


(pb/GeV) 


(%) 


(%) 




(%) 




30-40 


37.0 


1.47x 10~ x 


8 


+15/- 


11 +17/- 


14 


40-50 


44.8 


3.06x10 


4 


+14/- 


10 +15/- 


11 


50-60 


54.5 


1 A A . 1 A 1 

1.44x10 


4 


ill/ 

+11/- 


-9 


+ 12/- 


10 


art 7n 


04.0 


t no w 1 n — 2 


r 



111/ 

+11/- 


-9 


+12/- 


11 


70-80 


74.6 


4.21xl0~ 2 


7 


+ 14/- 


12 +16/- 


14 


80-90 


84.6 


2.57xl0 -2 


7 


+13/- 


11 +14/- 


12 


90-100 


94.8 


1.53xl0~ 2 


9 


+ 14/- 


13 +16/- 


15 


100-125 


110.9 


7.97xl0" 3 


6 


+ 12/- 


10 +14/- 


12 


125-150 


136.2 


2.88xl0~ 3 


7 


+ 15/- 


14 +16/- 


16 


150-200 


170.4 


1.27xl0 -3 


7 


+15/- 


13 +16/- 


15 


200-350 


249.2 


2.66xl0~ 4 


8 


+15/- 


14 +17/- 


17 


350-500 


403.0 


2.70xl0 -5 


22 


+47/- 


47 +53/- 


52 




pV 


(pV) 


da Idpip 


5atat 








(GeV) 


(GeV) 


(pb/GeV) 


(%) 


(%) 




(%) 




0.0-2.5 


1.4 


6.27X10- 1 


4 


+9/- 


9 


+10/- 


10 


2.5-5.0 


3.6 


5.38xl0 _1 


12 


+9/- 


9 


+15/- 


15 


5.0-7.5 


6.2 


3.50xl0 _1 


14 


+10/- 


-9 


+17/- 


17 


7.5-10 


8.8 


3.47xl0 _1 


15 


+15/- 


10 


+22/- 


19 


10-12.5 


11.2 


2.35X10" 1 


12 


+12/- 


11 


+17/- 


17 


12.5-15 


13.7 


1. 77x10^ 


16 


+ 12/- 


11 


+20/- 


20 


15-20 


17.3 


1.26xl0 _1 


10 


+ 12/- 


11 


+16/- 


15 


20-25 


22.4 


6.99xl0~ 2 


8 


+ 12/- 


11 


+15/- 


14 


25-30 


27.4 


5.29xl0" 2 


10 


+ 12/- 


10 


+16/- 


15 


30-40 


34.8 


6.32xl0~ 2 


8 


+ 12/- 


11 


+14/- 


14 


40-50 


44.5 


5.04xl0" 2 


9 


+ 13/- 


13 


+16/- 


16 


50-60 


54.7 


2.53xl0~ 2 


13 


+ 13/- 


12 


+19/- 


19 


60-80 


67.9 


1.04xl0" 2 


12 


+12/- 


11 


+17/- 


17 


80-100 


87.7 


3.45xl0~ 3 


17 


+20/- 


20 


+26/- 


26 


100-120 


108.4 


1.19xl0" 3 


19 


+20/- 


19 


+28/- 


28 


120-170 


139.6 


4.75xl0" 4 


20 


+20/- 


20 


+29/- 


28 



and acceptance calculations and by the statistical uncer- 
tainty in data, taking into account statistical correlations 
with adjacent bins. The latter are estimated using an in- 
verted smearing matrix, following a procedure described 
in Ref. |4lj . The smearing matrix represents the detec- 
tor resolution function and relates each bin at the par- 
ticle level to the bins at the reconstruction level. It is 
constructed for each variable using the DDP MC events 
simulated with SHERPA. 

Figures [T]-[3] show a comparison of the measured differ- 
ential cross sections to the theoretical predictions from 

DIPHOX, RESBOS, 27NNLO, and SHERPA. The RESBOS 

predictions are valid only for the phase space limited 
by 9 < M 11 < 350 GeV. We take this into account 
in our calculations and compare RESBOS predictions to 
Af 77 measurements up to «250 GeV (see Table II), the 
last mass value below 350 GeV where the cross section 
is measured. Systematic uncertainties across the bins in 
the measured cross sections are largely (> 90%) corre- 
lated. A common normalization uncertainty of 7.4% re- 
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FIG. 1: (Color online) The differential cross section as a function of (a) M 77 , (b) p~jp , (c) A0 77 , and (d) |cos(9*| for the full 
A</> 77 region from data (black points) and theory predictions (curves) are shown in the upper plots. The lower plots show the 
ratio of data and diphox, resbos, and 27NNLO predictions to the SHERPA predictions. The inner line for the error bars in data 
points shows the statistical uncertainty, while the outer line shows the total (statistical and systematic added in quadrature) 
uncertainty after subtracting the 7.4% normalization uncertainty. 



suiting from luminosity and diphoton selection efficiency 
is not shown in the plots. The predictions from SHERPA, 
diphox and RESBOS arc computed using the CTEQ6.6M 
NLO PDFs l42ll and from 27NNLO using MSTW2008 
NNLO PDFs [4J]. The PDF uncertainty is estimated 
using diphox and the 44 eigenvectors provided with the 
CTEQ6.6M PDF set. They are found to be within (3- 
7)%. The renormalization /x^, factorization /j,f, and frag- 
mentation [if scales are set to (1r = fip = fit = M 77 . 
The uncertainty due to the scale choice is estimated us- 
ing diphox via a simultaneous variation by a factor of 
two of all scales relative to the default values and found 
to be about 10% for dc/c?M 77 and dcr/d\ cos 9* |, and a 
maximum of (20-28)% for dajdpip at high pj 7 and for 



dcr I dA.(j>y~/ at low A0 77 . All theoretical predictions are 
obtained using diphoton event selection criteria equiva- 
lent to those applied in the experimental analysis (as are 
those used for the acceptance calculation). In particular, 
the photon is required to be isolated by < 2.5 GeV. 
For diphox, RESBOS, and 27NNLO, p^* is computed at 
the parton level. The cross sections from diphox, RES- 
BOS and 27NNLO are corrected for effects stemming from 
multiple parton interactions and hadronization, while for 
SHERPA these effects arc handled within the software 
package. These corrections are estimated using diphoton 
events simulated by pythia with Tunes A and SO [37I ]. 
The corrections vary within (4-6)% as a function of the 
measured kinematic variables and are consistent for both 
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M n (GeV) tt p; T (GeV) * |cose*| 



FIG. 2: (Color online) The differential cross section as a function of (a) M 77 , (b) p7p , and (c) |cos#*| for the A0 77 < n/2 
region. The notations for points, lines and shaded regions are the same as in Fig. [T] 




M n (GeV) tt p; T (GeV) tt |cose*| 



FIG. 3: (Color online) The differential cross section as a function of (a) M 77 , (b) pP , and (c) |cos#*| for the A0 77 > 7r/2 
region. The notations for points, lines and shaded regions are the same as in Fig. [1] 



tunes within 1%. 

Tables HTHVl show that the cross sections in the A0 77 > 
7r/2 region constitute, on average, about (85-90)% of 
the cross sections for the full A(/> 77 range. From the 
sub-tables for the pp variable, we observe that at pP< 
25 GeV, the cross sections are fully dominated by the 
A0 77 > 7r/2 region, while starting from pP> 30 GeV, 
they are significantly dominated (by a factor of 2-4) by 
the A</> 77 < 7r/2 region. The shoulder-like structure ob- 
served in the pP distribution around 30 — 40 GeV should 
be mainly caused by the fragmentation photons coming 
from the A0 77 < 7r/2 region, and partially by higher- 
order (NLO and beyond) corrections (l9j . 

In general, none of the theoretical models considered 
here provides a consistent description of the experimental 
results in all kinematic regions. The SHERPA predictions 



are able to describe most of the phase space relatively well 
except for the low DDP mass region, very low A</> 77 , and 
with some tension in the | cos 8* | spectrum. A noticeable 
discrepancy between RESBOS and DIPHOX in some regions 
of the phase space is due to the absence of all-order soft- 
gluon resummation (pP close to zero and A0 77 close to 
it) and the fact that the gg — > 77 contribution is cal- 
culated only at LO in diphox (small M 77 ). However, 
RESBOS fails to describe M 77 , pP , and |cos#*| spectra 
in the A0 77 < ir/2 region, where the contributions from 
the fragmentation diagrams and higher-order corrections 
arc important. The processes with a parton-to-diphoton 
fragmentation taking place at low masses (M 77 < pP) 
are not included yet in any existing calculation The 
regions of phase space with a significant contribution 
from fragmentation photons (very low A</> 77 ) require ex- 
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TABLE III: The measured differential cross sections in bins 
of A</> 77 and | cos 9* \ . The columns <5 stat , <5 sy st and <5 tot rep- 
resent the statistical, systematic, and total uncertainties, re- 
spectively. 



A0 77 


(A<£ 77 ) 


da 1 dA(f>-yy 


<5stat 


(5syst 


Stot 


(rad) 


(rad) 


(pb/rad) 


(%) 


(%) 


(%) 




0.00-0.31 


0.17 


2.28 


22 


+ 12/- 


12 +25/- 


25 


0.31-0.63 


0.46 


1.93 


16 


+14/- 


13 +21/- 


21 


0.63-0.94 


0.79 


5.66xl0 _1 


12 


+21/- 


21 +25/- 


24 


0.94-1.26 


1.11 


4.09xl0 -1 


13 


+21/- 


19 +25/- 


23 


1.26-1.57 


1.42 


5.62xl0 _1 


20 


+ 18/- 


17 +27/- 


26 


1.57-1.88 


1.73 


6. 82x10^ 


11 


+ 16/- 


14 +20/- 


18 


1.88-2.20 


2.05 


1.04 


8 


+ 14/- 


13 +17/- 


15 


2.20-2.51 


2.37 


1.65 


11 


+ 14/- 


12 +18/- 


17 


2.51-2.67 


2.60 


3.57 


13 


+21/- 


11 +25/- 


17 


2.67-2.83 


2.75 


4.98 


7 


+ 13/- 


11 +14/- 


13 


2.83-2.98 


2.91 


1.08X10 1 


6 


+13/- 


-9 +15/- 


11 


2.98-3.14 


3.08 


2.75 xlO 1 


3 


+9/- 


8 +9/- 


9 




I cos e* | 


(| cos dcr/d|cos6»* 


5atat 


^syst 


Stot 






(Pb) 


(%) 


(%) 


(%) 




0.0-0.1 


0.05 


2.58X10 1 


6 


+9/- 


8 +11/- 


10 


0.1-0.2 


0.15 


2.22X10 1 


4 


+9/- 


9 +10/- 


10 


0.2-0.3 


0.25 


1.91X10 1 


5 


+10/- 


-9 +11/- 


10 


0.3-0.4 


0.35 


1.49X10 1 


5 


+9/- 


9 +11/- 


10 


0.4-0.5 


0.45 


9.91 


7 


+10/- 


-9 +12/- 


12 


0.5-0.6 


0.54 


5.20 


9 


+11/- 


10 +14/- 


14 


0.6-0.7 


0.64 


1.73 


12 


+17/- 


17 +21/- 


21 



tensive tuning of all of the considered event generators. 

In summary, we have presented measurements of dif- 
ferential cross sections of photon pair production in pp 
collisions at ^/s = 1.96 TeV as functions of M 77 , pip , 
A0 77 , and |cos#*| for photons with p T > 18(17) GeV 
and |ry| < 0.9 in the full A</> 77 range and for A0 77 < ir/2, 
A0 77 > 7r/2 separately. The cross sections are compared 
to the predictions made by the DIPHOX, resbos, 27NNLO 
and SHERPA MC generators. Overall, SHERPA provides 
the best description of the measured cross sections. The 
experimental results show discrepancies with all theoret- 
ical predictions in the regions of small A(/> 77 and small 
diphoton mass for A0 77 > n/2, with minor differences 
in the shapes of the |cos#*| distribution. The results 
are important for understanding of DDP production and 
tuning of modern generators to study SM phenomena 
and search for beyond the SM processes. 
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and FUNDUNESP (Brazil); DAE and DST (India); Col- 
ciencias (Colombia); CONACyT (Mexico); NRF (Ko- 
rea); FOM (The Netherlands); STFC and the Royal So- 
ciety (United Kingdom); MSMT and GACR (Czech Re- 
public); BMBF and DFG (Germany); SFI (Ireland); The 
Swedish Research Council (Sweden); and CAS and CNSF 



TABLE IV: The measured differential cross sections in bins 
of M 77 , pj 7 , and | cos#*| for A<^ 77 < n/2. The columns <5 a tat, 
Jayat and <5tot represent the statistical, systematic, and total 
uncertainties, respectively. 



M 77 


(M 77 ) 


da/dM 77 


<5atat 






<5tot 


(GeV) 


(GeV) 


(pb/GeV) 


(%) 


(%) 




(%) 




30-40 


34.3 


1.64x10"^ 


14 


+ 14/- 


14 +20/- 


20 


40-50 


44.8 


8.92X10 -3 


28 


+ 15/- 


14 +31/- 


31 


50-60 


54.6 


2.25xl0" 3 


23 


+23/- 


23 


+33/- 


33 


60-70 


64.6 


1.22xl0" 3 


41 


+25/- 


27 +48/- 


49 


70-90 


78.7 


5.60xl0~ 4 


30 


+ 14/- 


14 +33/- 


33 


90-200 


111.4 


5.44xl0~ 5 


41 


+30/- 


30 +51/- 


51 




77 
Pt 


(pV) 


do /dp T y 


Sst&t 


<^syst 




<5tot 


(GeV) 


(GeV) 


(pb/GeV) 


(%) 


(%) 




(%) 




25-30 


28.3 


5.89X10 -3 


44 


+30/- 


28 


+54/- 


53 


30-40 


35.8 


3.56xl0" 2 


23 


+14/- 


14 


+27/- 


27 


40-50 


44.5 


4.39xl0~ 2 


15 


+17/- 


17 


+22/- 


22 


50-60 


54.8 


1.72xl0~ 2 


18 


+14/- 


14 


+23/- 


23 


60-80 


67.8 


7.74xl0~ 3 


17 


+ 12/- 


12 


+21/- 


21 


80-100 


87.5 


2.70xl0 -3 


22 


+17/- 


17 


+28/- 


28 


100-120 


108.3 


7.07xl0~ 4 


17 


+22/- 


22 


+28/- 


28 


120-170 


140.5 


3.84xl0~ 4 


25 


+26/- 


26 +36/- 


36 



| cos 6* | 


{Icosfl" 


|) do/d\cosG*\ 


<5atat 


<^syst 


Stot 






(Pb) 


(%) 


(%) 


(%) 




0.0-0.1 


0.05 


2.64 


23 


+23/- 


23 +33/- 


33 


0.1-0.2 


0.15 


3.22 


13 


+18/- 


18 +23/- 


23 


0.2-0.3 


0.25 


3.71 


15 


+14/- 


14 +21/- 


21 


0.3-0.4 


0.34 


2.17 


17 


+ 12/- 


12 +21/- 


21 


0.4-0.5 


0.45 


1.09 


17 


+15/- 


15 +23/- 


23 


0.5-0.6 


0.54 


6.12X10" 1 


39 


+23/- 


23 +45/- 


45 


0.6-0.7 


0.63 


3.33X10 -1 


39 


+27/- 


27 +48/- 


48 



(China). 
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TABLE V: The measured differential cross sections in bins of 
M 77 , p"Jp , and |cos#*| for A<^ 77 > tt/2. The columns <5 s tat, 
(5 sys t and <5 to t represent the statistical, systematic and total 
uncertainties, respectively. 



M 77 


(M 77 ) 


da / dM 11 


Sstat 


<5syst 




Stot 


(GeV) 


(GeV) 


(pb/GeV) 


(%) 


(%) 




(%) 




30-40 


37.5 


1.31xl0 _1 


9 


+11/- 


-9 


+14/- 


13 


40-50 


44.8 


2.96xl0 _1 


5 


+9/- 


8 


+10/- 


-9 


50-60 


54.5 


1.43xl0 _1 


4 


+8/- 


8 


+10/- 


10 


60-70 


64.5 


8.06xl0" 2 


5 


+9/- 


9 


+10/- 


10 


70-80 


74.6 


4.10xl0~ 2 


8 


+9/- 


9 


+12/- 


12 


80-90 


84.6 


2.53xl0 -2 


7 


+12/- 


12 


+13/- 


13 


90-100 


94.7 


1.53xl0~ 2 


9 


+12/- 


12 


+14/- 


14 


100-125 


110.9 


7.97xl0" 3 


7 


+10/- 


10 


+12/- 


13 


125-150 


136.2 


2.82xl0" 3 


8 


+15/- 


15 


+17/- 


17 


150-200 


170.4 


1.26xl0" 3 


7 


+13/- 


13 


+15/- 


15 


200-350 


249.2 


2.65xl0" 4 


9 


+15/- 


15 


+17/- 


17 


350-500 


403.0 


2.70xl0 -5 


28 


+47/- 


47 


+55/- 


55 




Pt 


(pV) 


da/dp'^ 1 


<5stat 


<^syst 




Stot 


(GeV) 


(GeV) 


(pb/GeV) 


(%) 


(%) 




(%) 




0.0-2.5 


1.4 


6.34xl0 _i 


5 


+9/- 


9 


+10/- 


10 


2.5-5.0 


3.6 


5.42xl0 _1 


16 


+8/- 


8 


+18/- 


18 


5.0-7.5 


6.2 


3.52X10" 1 


20 


+9/- 


9 


+22/- 


22 


7.5-10 


8.8 


3.48xl0 _1 


14 


+12/- 


-9 


+19/- 


17 


10-12.5 


11.2 


2.35X10" 1 


17 


+11/- 


11 


+20/- 


20 


12.5-15 


13.7 


1.77X10" 1 


16 


+11/- 


11 


+20/- 


20 


15-20 


17.3 


1.26X10" 1 


11 


+10/- 


10 


+15/- 


15 


20-25 


22.4 


6.96xl0" 2 


9 


+10/- 


10 


+14/- 


14 


25-30 


27.3 


4.82xl0" 2 


11 


+12/- 


12 


+16/- 


17 


30-40 


34.2 


3.03xl0" 2 


8 


+10/- 


10 


+13/- 


13 


40-50 


44.4 


1.21xl0~ 2 


10 


+9/- 


9 


+14/- 


14 


50-60 


54.5 


5.93xl0" 3 


14 


+20/- 


20 


+24/- 


24 


60-80 


68.2 


2.18xl0~ 3 


14 


+16/- 


16 


+21/- 


21 


80-100 


88.3 


5.83xl0" 4 


2r> 


+29/- 


29 


+38/- 


38 


100-120 


108.5 


4.91xl0~ 4 


27 


+24/- 


24 


+36/- 


36 


120-170 


137.4 


1.13xl0 -4 


30 


+31/- 


31 


+43/- 


43 



cos <9*| 


(|cos0 


*|) dcr/d|cos6>* 


<5stat 


•^syst 




<5tot 






(Pb) 


(%) 


(%) 




(%) 




0.0-0.1 


0.05 


2.24X10 1 


6 


+8/- 


8 


+10/- 


10 


0.1-0.2 


0.15 


1.86X10 1 


5 


+8/- 


8 


+9/- 


9 


0.2-0.3 


0.25 


1.55X10 1 


6 


+9/- 


8 


+11/- 


11 


0.3-0.4 


0.35 


1.24X10 1 


6 


+8/- 


8 


+10/- 


10 


0.4-0.5 


0.45 


8.38 


8 


+9/- 


9 


+ 12/- 


12 


0.5-0.6 


0.54 


4.43 


10 


+11/- 


10 +15/- 


14 


0.6-0.7 


0.64 


1.04 


17 


+22/- 


22 +28/- 


28 



[11] M. Begel, Ph.D. Thesis, University of Rochester (1999), 

FERMILAB-THESIS-1999-05. 
[12] C. Albajar et al. (UA1 Collaboration), Phys. Lett. B 209, 

385 (1988). 

[13] J. Alitti et al. (UA2 Collaboration), Phys. Lett. B 288, 

386 (1992). 

[14] V. M. Abazov et al. (DO Collaboration), Phys. Lett. B 

690, 108 (2010). 
[15] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. Lett. 

107, 102003 (2011). 
[16] G. Aad et al. (ATLAS Collaboration), Phys. Rev. D 85, 



012003 (2012). 

[17] S. Chartrchyan et al. (CMS Collaboration), J. High En- 
ergy Phys. 1201, 133 (2012). 

[18] T. Binoth, J. Ph. Guillet, E. Pilon, and M. Werlen, Eur. 
Phys. J. C 16, 311 (2000). 

[19] T. Binoth, J. Ph. Guillet, E. Pilon, and M. Werlen, Phys. 
Rev. D 63, 114016 (2001). 

[20] J. C. Collins and D. E. Soper, Phys. Rev. D 16, 2219 
(1977). 

[21] The polar angle 6 is defined with respect to the positive z 
axis, which is along the proton beam direction, while the 
azimuthal angle <j> is defined with respect to the x axis. 
Pseudo-rapidity is defined as r\ = — ln[tan(6'/2)]. Also, 
|?7det| and |</>det| are the pseudorapidity and the azimuthal 
angle defined with respect to the geometric center of the 
DO detector. 

[22] P. Nadolsky, C. Balazs, E. Berger, C.-P. Yuan, Phys. Rev. 

D 76, 013008 (2007). 
[23] C. Balazs, E. Berger, S. Mrenna, and C.-P. Yuan, Phys. 

Rev. D 57, 6934 (1998). 
[24] S. Catani, L. Cieri, D. Florian, G. Ferrera, and M. Grazz- 

ini, Phys. Rev. Lett. 108, 072001 (2012). 
[25] T. Gleisberg et al., J. High Energy Phys. 02, 007 (2009). 

We use sherpa v.1.2.2. 
[26] S. Hoche, S. Schumann, and F. Siegert, Phys. Rev. D 81, 

034026 (2010). As suggested by the sherpa authors, we 

choose Qcut = 10 GeV as the ME-PS matching scale and 

parameter D = 0.4 taken to be of the size of the photon 

isolation cone. 

[27] S. Catani, B. R. Webber, and G. Marchesini, Nucl. Phys. 

B 349, 635 (1991). 
[28] A. Banfi, G. P. Salam, and G. Zanderighi, J. High Energy 

Phys. 05, 073 (2005). 
[29] S. Schumann and F. Siegert, private communication. 
[30] S. Catani and M. Grazzini, Phys. Rev. Lett. 98, 222002 

(2007). 

[31] V. M. Abazov et al. (DO Collaboration), Nucl. Instrum. 

Methods Phys. Res. A 565, 463 (2006). 
[32] M. Abolins et al, Nucl. Instrum. Methods in Phys. Res. 

Sect. A 584, 75 (2008). 
[33] R. Angstadt et al, Nucl. Instrum. Methods in Phys. Res. 

Sect. A 622, 298 (2010). 
[34] R. Brun and F. Carminati, CERN Program Library Long 

Writeup W5013, 1993 (unpublished). 
[35] T. Andeen et al, FERMILAB-TM-2365 (2007). 
[36] V. M. Abazov et al. (DO Collaboration), Phys. Rev. Lett. 

102, 231801 (2009). 
[37] T. Sjostrand et al., Comput. Phys. Commun. 135, 238 

(2001) . 

[38] V. M. Abazov et al. (DO Collaboration), Phys. Lett. B 

653, 378 (2007). 
[39] R. Hamberg, W. L. van Neerven, and T. Matsuura, 

Nucl. Phys. B359, 343 (1991) [Erratum-ibid. B644, 403 

(2002) ]. 

[40] C. Buttar et al. , larXiv:0803.0678l [hep-ph], section 9. 
[41] V. M. Abazov et al. (DO Collaboration), Phys. Rev. D 

83, 052008 (2011). 
[42] J. Pumplin et al, J. High Energy Phys. 0207, 012 (2002) 

and D. Stump et al, J. High Energy Phys. 0310, 046 

(2003) . 

[43] A. D. Martin, W. J. Stirling, R. S. Thorne, and G. Watt, 
Eur. Phys. J. C 63, 189 (2009). 



